J/iIj suppression in p— A collisions from parton energy loss in cold QCD matter 
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The effects of energy loss in cold nuclear matter on J/^p suppression in p-A collisions are studied. 
A simple model based on first principles and depending on a single free parameter is able to reproduce 
J/ip suppression data at large and at various center-of-mass energies. These results strongly 
support energy loss as a dominant effect in quarkonium suppression. They also give some hint 
on its hadroproduction mechanism suggesting color neutralization to happen on long time-scales. 
Predictions for J/tp and T suppression in p-Pb collisions at the LHC are made. 
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The spectacular results on jet production in Pb-Pb 
collisions at the LHC (see e.g. [H, 01) find a natural ex- 
planation in terms of parton energy loss in quark-gluon 
plasma. Despite the wealth of data accumulated so far 
at RHIC and LHC, the in-depth understanding of en- 
ergy loss processes remains far from complete (see Q 
for a discussion). Recently, new scaling properties have 
been identified for the induced gluon radiation spectrum, 
dl/dw, of hard processes where a color charge undergoes 
small angle scattering through a static medium (nuclear 
matter or quark-gluon plasma) 4] . 

The goal of this Letter is to explore the phenomeno- 
logical consequences of these results on J/ ?/> as well as T 
suppression, 
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in p-A and tt-A collisions (we keep in the following the 
notations "J/i/j" and "p-A" for clarity). This observable 
allows for probing energy loss in nuclear matter, which 
is a well-controlled medium as opposed to an expanding 
quark-gluon plasma which dynamics is more complex. 
Moreover, understanding quarkonium suppression in p- 
A collisions is a prerequisite in order to interpret quanti- 
tatively the measurements performed in heavy-ion colli- 
sions at the LHC 0-0] • It is striking that there is no con- 
sensus yet on the origin of the significant J/tp suppression 
reported at large rapidity in p-A collisions, from SPS to 
RHIC (8l-[lo|. despite many theoretical attempts (see 11 1 
for a review). 

As is well-known, quarkonium hadroproduction in el- 
ementary p-p collisions is not fully understood. In or- 
der to study quarkonium nuclear suppression in the most 
model-independent way, the p-p quarkonium production 
cross section will be taken from experiment. We will only 
assume that the heavy-quark QQ pair is produced in a 
compact color octet state, within the hard process time- 
scale ihj and remains color octet for a time much longer 



than th- In quarkonium production models where color 
neutralization is a soft, non-perturbative process, this as- 
sumption holds at any . In the Color Singlet Model, we 
expect this assumption to be founded at large enough , 
where the gluon emission required for color neutralization 
is constrained to be semi-hard (or even softish) by energy 
conservation. With this assumption, at sufficiently large 
quarkonium energy E in the target rest frame, quarko- 
nium hadroproduction looks like small angle scattering 
of a color charge [2^. The associated soft gluon radia- 
tion spectrum is thus similar to the (non-abelian) Bethe- 
Heitler spectrum of an asymptotic charge, and depends 
on the amount of transverse momentum kick q± to the 
charge. The typical q± is expected to be larger in p-A 
than in p-p collisions due to transverse momentum nu- 
clear broadening Aq^. As a result, the medium-induced 
radiation spectrum is similar to the Bethe-Heitler spec- 
trum with replaced by Ag^ (see Q), 
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In the following we shall take ag — 0.5, A — Aqqd = 
0.25 GeV, — 1 GeV in the transverse mass Mj_ = 
+ pI , and A/ = 3 GeV (M = 9 GeV) for the mass 
of a compact cc (65) pair. 

This leads to an average medium-induced radiative loss 
scaling as the quarkonium energy, AE oc E. In the limit 
A^ < Aq]_ < Ml we have 

fE 



AE= I dLooj—^N.asy^^E. 

duj Mj_ 



(3) 



The scalingAiJ cx E was first postulated in [l^l (also re- 
visited in [13|) yet this assumption was not motivated and 
the parametric dependence on L and Af arbitrary (and 
different from ([3])). In Ref. jl4l |. a bound on medium- 
induced energy loss was derived, AE < E^ , but in a 
specific setup where the nuclear broadening of the final 
tagged particle was neglected. 
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The starting point of the model consists in expressing 
the J/V" differential production cross section dcr/da;^ in 
p-A collisions simply as that in p-p collisions, with a 
shift in Xp accounting for the energy loss e incurred by 
the octet cc pair propagating through the nucleus, 



1 day 

A dx^ 



{x,)= I AeV{e)^^{x,+5x,{e)), 

(4) 

where x^ is defined as the longitudinal momentum frac- 
tion between the J /ip and projectile proton in the cm. 
frame of an elementary p-N collision (of energy \fs). In 
the limit ^/s 3> mp (with ni-p the proton mass), it reads 



' E s 



(5) 



where E and E^ ~ s/(2mp) are respectively the cc pair 
and projectile proton energies in the nucleus rest frame. 
We now describe the various ingredients in Eq. (jH): 
(i) The differential p-p cross section is determined from 
a fit of p-p data and can be conveniently parameterized 
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pp 
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(x^) oc {l~xT/x' ; x' = ^xl + 4Mi/s . (6) 



The exponent n is extracted from p-p data taken at the 
same cm. energy (whenever possible) as the p-A mea- 
surements discussed in this Letter (the normalization pa- 
rameter being irrelevant for our purpose, see ([Ij); 
(ii) The shift Sx^, (e) is defined by 



Xp{E) + 5xp{e) = Xp{E+e) 
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where E is obtained directly from ([S]). Note that at large 
Xp ^ M±/y/s, we have E c^ix^E-p and ^^^(e) ~ s/E-p] 
(iii) The average over e in Eq. Q is performed using the 
energy loss probability distribution, or quenching weight 
V{e). As a simple and physically sound choice for a nor- 
malized distribution, we take |15| 
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The amount of medium-induced gluon radiation, and 
hence the strength of J /ip suppression in p-A collisions, 
is controlled by Ag^ in Eq. For a path length L 

travelled across the target (proton or nucleus), it reads 



Ag5_(i) = QA L - qp Lp 



(9) 



where the average path length is given by i = | ro A^/"^ 
(ro = 1.12 fm), assuming the hard process to occur uni- 
formly in the nuclear volume. In 0, qa (resp. qp) stands 
for the transport coefficient in the nucleus A (resp. in the 



proton). It is related to the gluon distribution G{x) in a 
target nucleon as (l6| 



q(x) 



m ~ 1 



pxG{x, qL) , 
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where p is the target nuclear density. Apart from the scal- 
ing violations in the running of as and in the evolution 
of the gluon density, both neglected since qL <1 GeV^, 
the L-dependence of q mainly enters via the typical x 
value at which xG{x) should be evaluated. When the 
hard production time ih <SC L, it is estimated to be 
X ~ (2r7ipL)~^ = xo dl]- When th 3> L, the hard sub- 
process is coherent over the whole nucleus, and in this 
case we expect x ^ X2 [l5| , where X2 is the target parton 
momentum fraction, X2 = {—x^ + x')/2 when assum- 
ing a 2 — > 1 subprocess kinematics. Using the power- 



law behavior xG{x 
data 



-0.3 



suggested by fits to HERA 



17[ , (7 is thus given by [26 
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; X = min(a::o,X2) . 
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In (fTTj) . the transport coefficient qo = q{x ~ 10^^) is the 
only free parameter of the model. 

Besides energy loss effects, other mechanisms might af- 
fect J/ip suppression in nuclei. At small energy E, when 
the J/ip hadronization time t^ = jt^ < L, J/ip^s are 
produced inside the nuclear medium, and consequently 
might suffer inelastic interaction with the target nucleus, 
or nuclear absorption. In the figures below, we indicate 
by an arrow the typical value of below which this 
starts to happen, assuming ~ 0.3 fm. Another effect 
is the expected saturation of the nuclear gluon density 
at small X2, leading to an additional J/^ suppression in 
high-energy p-A collisions. The associated suppression 
is a scaling function of the saturation scale Qs, which can 
be simply parameterized as [l8| 



{h + Ql{x2.L)r- 



(12) 



In order to make reliable predictions at RHIC and LHC, 
where saturation effects might be important, the J /ip nu- 
clear production ratio will be determined assuming en- 
ergy loss effects, R^'a^^ from Eq. (jl]), with and without 
saturation effects, 

(i) RpA = ^pA°'''' : 

(ii) RpA - X sT/sf . 



The saturation scale appearing in (|12|) is determined con- 
sistently in our model through the relationship jl9| 



Ql{x,L)^q{x)L, 



(13) 



where q{x) is given by (jlll) . The inclusion of saturation 
effects thus does not require any additional parameter. 
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Figure 1: E866 J/^ suppression data [1] in p-W (left) and 
p-Fe (right) collisions compared to the energy loss model. 



The only parameter of the model, the transport co- 
efficient (Jq, is determined by fitting the J/ip suppres- 
sion measured by E866 Q in p-W over p-Be collisions 
at = 38.7 GeV in the [0.2-0.8] x^-range f27|. This 
choice is motivated by the fact that these data are the 
most precise measurements performed so far and cover a 
wide range in Xp. The fit gives qo — 0.09 GeV^/fm as- 
suming energy loss effects only, and qo — 0.05 GeV^/im 
when saturation effects are also included. The result of 
the fit in these two cases is shown in Fig. [1] (left panel); 
the agreement is excellent in the fit range, while a slight 
disagreement is observed below x^, ^0.1, where nuclear 
absorption is expected to play a role. The successful de- 
scription of J/ ip suppression in iron, Rpc/Bc (right panel) , 
at the same energy fully supports the atomic mass depen- 
dence of the model. Note also that the values for the 
transport coefficient, qq — 0.05-0.09 GeV^/fm, would 
correspond to Ql{x = IQ-^) = 0.08-0.15 GeV^ using 
([T3| . which is consistent with (yet slightly smaller than) 
estimates based from fits to F2 DIS data (20| . 

Data taken at lower ^/s or smaller also com- 

pared to the model. As can be seen in Fig. [2] the agree- 
ment is very good, both in shape and magnitude, over a 
very wide range in x^ . We also remark that as expected, 
the model tends to slightly underestimate the suppression 
in the x^ domain where J/ip^s might hadronize inside the 
nucleus and suffer nuclear absorption; see in par ticular 
the comparison with NA60 [211] and HERA-B 0| data. 
It is also remarkable that the model is able to reproduce 
the different magnitude of suppression in p-A and 7r~- 
A collisions reported by NA3 [8| (Fig. ^ upper panels). 
This difference cannot be understood within nuclear ab- 
sorption models (since factorization is usually assumed 
between the hard production process and final state inter- 
action) nor can it be explained by nuclear PDF (nPDF) 
effects, unless the nPDF to proton PDF ratios for valence 
quarks and for gluons, probed respectively in 7r~-A and 
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Figure 2: NA3 [1], NA60 [21| and HERA-B ^ J/^ suppres- 
sion data in p-A (7r~-A) collisions compared to the energy 
loss model. 



p-A collisions, prove completely different. In our picture, 
the smaller J/ip suppression in tt~-A collisions naturally 
arises from the flatter differential cross section, n^p = 1.5 
vs. Upp = 4.3 (see (O) at NA3 energies, a feature which 
can be explained from the slope of the PDF in a pion and 
in a proton, respectively. For completeness, the predic- 
tions including saturation are also shown as dashed lines 
in Fig. [21 as expected saturation effects are tiny at these 
energies. The slight differences are essentially due to the 
smaller transport coefficient used (go = 0.05 GeV^/fm) 
when saturation is included. 
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Figure 3: PHENIX J/-(/j suppression data [iQl] in d-Au colli- 
sions compared to the energy loss model. 

The predictions in d-Au collisions at RHIC, ^/s = 



4 



200 GeV, are shown in Fig. |3] in comparison with 
PHENIX data Energy loss effects are able to re- 

produce J/ij} suppression at positive rapidities. However, 
some disagreement is observed in the negative y bins for 
which nuclear absorption is expected to play a role; it 
has also been conjectured that a depletion of the gluon 
nPDF at large X2 might explain the trend of the data at 
backward rapidities [23| . We also note that the disagree- 
ment is reduced when saturation effects are included. Fi- 
nally, arguing about the possible slight disagreement ob- 
served around mid-rapidity might be premature, given 
the present experimental uncertainties. 
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Figure 4: J/ip and T suppression expected in p-Pb collisions 
at the LHC. 

Last, the rapidity dependence of J/i/; suppression in p- 
Pb collisions at the LHC (taking -/s — 5 TeV) is shown in 
Fig. 01 Even though the inclusion of saturation effects is 
expected to yield a stronger J/ %() suppression, it is some- 
how compensated by the use of a smaller transport co- 
efficient; as a consequence predictions with (dashed line) 
and without (solid) saturation are actually rather simi- 
lar, except in the negative y-bins. As can be seen from 
the arrow in Fig.|31 J/i) hadronization should take place 
outside the nuclear medium above y > — 5; nuclear ab- 
sorption should thus play little or no role at the LHC. 
At forward rapidities, J/tj) suppression becomes rather 
large, e.g., RpPh — 0.7-0.8 at y = 1 down to i?ppb ^ 0.5 
at y > 4. Fig. m also shows the predicted T suppression as 
a dash-dotted line [2^ . Because of the mass dependence 
of energy loss, AE (x Mj_~^ , it is expected to be smaller 
than that of J ftp yet not negligible, e.g., R^pi, — 0.8 
at y = 3. These predictions can be compared to the 
future measurements by the ALICE and LHCb experi- 
ments during the p-Pb run scheduled in 2012. 

In summary, an energy loss model (supplemented by 
saturation effects) based on first principles has been pre- 
sented. Once the transport coefficient is determined from 
a subset of E866 data, it is able to reproduce nicely all 
existing J/ijj measurements in p-A collisions. In particu- 



lar the dependence of J/ ifj suppression on Xj, , for various 
atomic mass A and center-of-mass energies ^/s, is well 
accounted by the model. These results strongly support 
energy loss in cold nuclear matter to be the dominant 
effect of J/V' suppression, at least at large rapidities. It 
also qualitatively explains why no sup pre ssion in inelas- 
tic J/ijj electroproduction is observed 2Jj, since Eq. ([3]) 
does not apply to this case Q . Finally, the agreement be- 
tween the data and our model predictions supports the 
assumption of a long-lived color octet QQ pair. As a per- 
spective, we plan to investigate energy loss effects on J/ip 
suppression in heavy-ion collisions as well as to explore 
their consequences on other hard processes. 
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Lacking T p-p data, we assume n 



in Eq. (O. 



